Magnetic phase diagram of CuO 
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High resolution ultrasonic velocity measurements have been used to determine the temperature - 
magnetic- field phase diagram of the monoclinic multiferroic CuO. A new transition at Tn3 = 230 K, 
corresponding to an intermediate state between the antiferromagnetic non-collinear spiral phase 
observed below Tn2 = 229.3 K and the paramagnetic phase, is revealed. Anomalies associated with 
a first order transition to the commensurate collinear phase are also observed at Tjvi = 213 K. 
For fields with B || b, a spin- flop transition is detected between 11 T - 13 T at lower tempera- 
tures. Moreover, our analysis using a Landau- type free energy clearly reveals the necessity for an 
incommensurate collinear phase between the spiral and the paramagnetic phase. This model is also 
relevant to the phase diagrams of other monoclinic multiferroic systems. 

PACS numbers: 75.10.-b, 75.30.Kz, 75.85.+t, 75.30.Gw 



Multiferroic phenomena have been a subject of intense 
interest in recent decades arising from opportunities to 
explore new fundamental physics as well as possible tech- 
nological applications [1-3]. Coupling between different 
ferroic orders has been proven to be driven by several dif- 
ferent types of mechanisms. In particular, multiferroics 
with a spiral spin-order-induced ferroelectricity have re- 
vealed high spontaneous polarization and strong magne- 
toelectric coupling [4, 5]. Cupric oxide (CuO), the sub- 
ject of this letter, was characterized as a magnetoelectric 
multiferroic four years ago when it was shown that its fer- 
roelectric order is induced by the onset of a spiral antifer- 
romagnetic (AFM) order at an unusually high tempera- 
ture of 230 K [3]. Thus far, two AFM states have been re- 
ported, a low temperature (Tjvi ~ 213 K) AF1 commen- 
surate collinear state with the magnetic moments along 
the monoclinic b axis and an AF2 incommensurate spi- 
ral state with half of the magnetic moments in the ac 
plane (Xjv2 ^ 230 K) [3, 6, 7]. However, the authors of 
the neutron diffraction measurements [6] questioned the 
possibility of having a direct condensation from a param- 
agnetic (PM) phase to a spiral magnetic phase. Despite 
this remark, a recent Landau theory [8], as well as several 
Monte-Carlo simulations [9, 10], appear to support this 
sequence of magnetic orderings. 

Encouraged by recent experiments on other multifer- 
roic systems using ultrasonic measurements [11], we mea- 
sured the temperature and field dependence of the veloc- 
ity of transverse modes in order to determine the mag- 
netic phase diagram of CuO. A new transition is detected 
at T/V3 = 230 K just above the AF2 spiral phase ob- 
served at Tn2 = 229.3 K, while the first order transition 
is observed at T/vi = 213 K. Furthermore, dielectric con- 
stant measurements confirm that only the spiral phase 
(between T/yi and T/V2) supports a spontaneous electric 
polarization. In addition, we report on a spin-flop tran- 



sition in the low temperature AF1 collinear phase when 
B || b. Thus, based on these findings, a new magnetic- 
field vs temperature phase diagram is proposed for CuO. 

In order to elucidate the possible nature of the AFM 
states observed in CuO, a non-local Landau- type free en- 
ergy is also developed for CuO and similar monoclinic 
multiferroics. This approach has been very successful in 
explaining the magnetic phase diagrams of other multi- 
ferroic systems [12-14]. In contrast with the conclusions 
of Refs. [8-10], our analysis based on rigorous symmetry 
arguments indicates that there must be a collinear in- 
termediate phase (AF3) between the paramagnetic and 
spiral AF2 states. Such a phase has been shown, both 
theoretically and experimentally, to occur in other geo- 
metrically frustrated antiferromagnets where symmetry 
allows for uniaxial anisotropy at second order [14, 15]. 
Finally, we compare the model predictions to the B-T 
phase diagram of CuO obtained using ultrasonic velocity 
data. Similarities with other multiferroic systems such 
as MnW0 4 , AMSi 2 6 , RMn0 3 , RMn 2 5 , and Ni 3 V 2 8 
are also noted. 

For the purpose of this study, a CuO sample was grown 
using a floating zone technique as described in Ref. [3]. 
A single crystal was cut with faces perpendicular to the 
monoclinic axes a*, b* = b, and c* (4x4x3 mm 3 ). The 
sample was then polished to obtain parallel faces. For 
velocity measurements, plane acoustic waves were gen- 
erated using 30 MHz LiNb03 piezoelectric transducers 
bonded to opposite faces. Using an ultrasonic interferom- 
eter, which measures the phase shift and the amplitude 
of the first elastic transmitted pulse, high-resolution rela- 
tive velocity variations (AV/V ~ 1 ppm) were achieved. 
Experimental data presented here were all obtained us- 
ing the velocity of transverse waves V a * [c*] propagating 
along the a* axis and polarized along c*, with the mag- 
netic field applied along the easy magnetic axis of CuO 



2 




FIG. 1. Temperature dependence of the dielectric constant 
€b (in green) and the relative velocity variations of transverse 
mode V a *[c*] measured at different fields with B || b. 



(b axis). Simultaneous capacitance measurements were 
carried out using an AH 2550A Ultra Precision 1kHz Ca- 
pacitance Bridge to identify which of these phases are 
ferroelectric. For that purpose, electrodes were mounted 
on faces perpendicular to the b axis in order to determine 
the dielectric constant e b . 

Fig. 1 shows the temperature dependence of the rela- 
tive sound velocity variations (AV/V) for B || b. At zero 
field, the anomaly observed at T/vi = 213 K (see inset of 
Fig. 1) coincides very well with the onset of a commen- 
surate collinear antiferromagnetic state. Our high reso- 
lution velocity measurements also reveal two anomalies 
at T N2 = 229.3 K and T N3 = 230.0 K near the stabiliza- 
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FIG. 2. Magnetic phase diagram of CuO for B|| b. Inset 
shows the relative velocity variation of V a *[c*] as a function 
of the field for T = 125 K. 



tion of a spiral order previously determined by neutron 
diffraction and susceptibility measurements [3, 6], which 
were thought to occur at a single transition. At higher 
fields, the amplitude of the step like variation observed 
at 229.3 K, as well as the temperature difference between 
T/V2 and T/v 3 increases, confirming the existence of a new 
intermediate magnetic order AF3. This finding is sup- 
ported by dielectric measurements also shown in Fig. 1. 
Notice that, as the stability range of the intermediate 
phase is small (AT ~ 0.7 K), velocity and dielectric data 
have been collected simultaneously to avoid any ambi- 
guity regarding the actual critical temperatures. Thus, 
as shown in Fig. 1 (for B = and 7 T), the anomaly 
observed on the dielectric constant e b coincides very well 
with T/V2 determined using velocity data, while no varia- 
tion is noticeable at T7V3. These results also indicate that 
the new phase AF3 is not ferroelectric, while magneto- 
electric coupling exists for the AF2 phase. 

We present in Fig. 2 the magnetic phase diagram of 
CuO determined up to 16 T using ultrasonic velocity 
measurements for B || b. The inset of Fig. 2 shows the 
field dependence of the velocity which displays a min- 
imum around 11 T for T = 125 K. As the magnetic 
moments are known to be parallel to the field in the 
AF1 commensurate collinear state [3, 6], we attribute 
this anomaly to a spin- flop transition [16]. In summary, 
while the critical temperatures Xjvi, Tiv 2 , and T/V3 are 
weakly field dependent, the spin- flop critical field Hsf 
increases with temperature. At 10 K, Hsf = 11 T and 
increases slowly up to 13.5 T at T/vi, in good agreement 
with magnetic susceptibility measurements performed on 
powder samples [17]. 

Since no neutron scattering data exists for the HF1 and 
AF3 states, we develop a Landau- type model in order to 
elucidate the nature of these new magnetic orders [14, 
15]. The integral form of the free energy is expanded in 
powers of the nonlocal spin density s(r) defined in terms 
of a uniform field-induced magnetization m and a spin 
polarization vector S modulated by a single wave vector 
Q describing the long-range magnetic order (Eq. (6) of 
Ref. [14]). Within the present model, the value of Q 
can be determined by simply considering the isotropic 
quadratic contribution 

F21 = ^2 J dridr 2 A(ri - r 2 )s(ri) • s(r 2 ), (1) 

which leads to F 2 i = \Am 2 + A Q S 2 where A Q = aT + 
Jq , with Jq being the Fourier transform of the exchange 
integral J(R). Considering the C-type monoclinic cell 
with four Cu 2+ magnetic ions, we obtain 

J(Q) = 2 [JiA(Q) + J 2 / 2 (Q) + J 3 / 3 (Q) + J 4 / 4 (Q)] 
/i(Q) = cos(7n? a - 7rq c ) 

/2(Q) = cos (7rq a + 7r<7 c ) (2) 
/3(Q) = cos (7rq a - 7rq b ) + cos (7rq a + 7rq b ) 
/ 4 (Q) = cos (7rq b - 7rq c ) + cos (7rq b + 7n? c ), 



where J\ and J2 represent the nearest-neighbors (NN) ex- 
change interactions along the AFM-chain (sites 2-3) and 
the coupling between chains (sites 1-4) on the same plane 
normal to b, respectively, and J3 and J4 represent the ex- 
change interactions along a (sites 1-2) and c (sites 1-3) 
between ions on different planes (see Fig. 4). The value of 
Q is then obtained by finding the extrema of Jq (Eq. 2) 
as a function of the exchange interactions. Results of our 
numerical algorithm are summarized in the J2 — J3 phase 
diagram shown in Fig. 3 for AFM chains (Ji = 1). For 
different J4 values, we obtain three phases: an incom- 
mensurate phase with Qjcm — [#a,0,g c ] (left side) and 
two commensurate phases (top and bottom right side). 
Depending on the sign of J3 relative to J4, the commensu- 
rate wave vector is either Q C m ~ [100] or Qcm ~ [001]. 
More interestingly, with J3 = J4 = we obtain the ex- 
pected commensurate wave vector Qcm = [\ — \) for 
J2 < (dash line in Fig. 3). Moreover, an ICM state 
with a modulation vector comparable to that of the ex- 
perimental value Qicm = [0.506 — 0.483] is stabilized 
whenever J3 and/or J4 are non-zero but small relative 
to Ji (for example, J2/J1 = —0.3, J3/J1 = 0.017, and 
J4/J1 = leading to Jq/ J\ = —2.6). These relative val- 
ues are also in good agreement with estimates obtained 
by density functional theory [9, 18, 19] and are consistent 
with the quasi- ID magnetic character of CuO. 

In addition to the usual isotropic second order ex- 
change term, we also consider anisotropic contributions. 
Considering the symmetry of monoclinic crystals (C2/c), 
we identified three invariants, written in single- ion form 
as 

F2A = ^j [ [D y (v)s y (r)s y (v) +D z (r) Sz (r) Sz (r) 
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FIG. 3. J3 — J 2 phase diagram for different values of J a with 
Ji = 1. One incommensurate phase with Q ICM = [qa,0,q c ] 
(left side) and two commensurate phases (Qcm = [100] and 
Qcm = [001]) are obtained. 
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FIG. 4. Spin configurations in a magnetic cell of 8 ions (red 
and orange circles). Red circles represent magnetic ions at b 
= 1/2. The +/- symbols represent spins in/out of the page. 
When no direction is specified (as in AF3 and HF3) , spins on 
these sites are not ordered. 

+ D xz (r)s x (r)s z (r)]dr . (3) 

While D y can be used to set the magnetic easy axis along 
b, the other terms are necessary in order to define the 
direction of the moments in the ac plane. Furthermore, 
to account for non-collinear spin configurations, we define 
S = Si + i S2, with 

51 = 5cos/3[cos7 y + sin 7 fa], (4) 

52 = £ sin /3[cos 6 pi + sin #(cos 7 y + sin 7 p 2 )] , 

where p\ and p2 are two orthogonal unit vectors nor- 
mal to the easy axis, y || b. Thus, the direction of the 
moments in the ac plane is accounted for by defining 
the unit vectors p\ and p2 relative to the lattice vectors, 
pi = cos a x + sin a z and p2 = —sin a x + cos a z. 
As shown in Fig. 4, the parameter a represents the an- 
gle between the ac plane component of S relative to the 
monoclinic axis a || x. After integration, all second-order 
contributions for m || H || y reduce to 

F total = l^ om 2 + AqS 2 _ l DyQm 2 _ D yQ \S y \ 2 

- D zQ \S z \ 2 + D xzQ S x S z - H m. (5) 

Adopting the same approach for the fourth-order 
isotropic term, we obtain 

F 4I = B X S A + i# 2 |S • S| 2 + i# 3 m 4 + 2£ 4 |m • S| 2 

+ B 3 m 2 S 2 + ^Bu[(S • S) 2 + c.c.]A 4 q,g • (6) 

Note the umklapp term A4q ; g, arising directly from the 
lattice periodicity [12]. This term is crucial in order to 
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account for the first order phase transition observed at 
Tni in CuO where a commensurate collinear state is sta- 
bilized. 

The free energy, F = F 2 i + F 2 a + Fu, with A Q = 
a(T — Tq) and A — D yQ = a(T — T ), is then numerically 
minimized. As in Ref. [15], most coefficients are set using 
analytical solutions associated with phase boundaries of 
second order transition. For example, setting Tq = 1.18, 
D yQ = 0.02, B x = 0.103, and B 2 = 0.011, reasonable 
values for the critical temperatures at zero field (T/V3 = 
1.2 and T/V2 = 1.12). We also set D z q = 0.01 as we must 
have D z q < D y Q, while the direction of the spins in the 
ac plane (a exp ~ 70°) [1] is used to determine the ratio 
D xz q/D z q = —0.42. The last coefficients are determined 
using the temperature of the multicritical point (where 
Tn2 and Tns boundaries meet) and the maximum field 
at T = K. From this exercise, we find B3 = 0.063 and 
B4 = 0.013 while B$ = 0.1 was set arbitrarily. Finally, 
Bu = 0.035 is used to obtained T N1 = 0.77. 

Fig. 5 shows the magnetic phase diagram obtained 
from minimization of the free energy. For comparison, 
we also present results obtained without the anisotropic 
terms D z q and D xz q (dotted lines). Depending on the 
scenario considered, we obtain 5 or 6 magnetic phases 
illustrated in Fig. 4, described by the order parameters 
listed in Table I. At zero field, both models (with and 
without D z and D xz ) predict the same phase sequence, 
consistent with our experimental observations shown in 
Fig. 2. At low temperatures, a collinear phase AF1 with 
the moments along b is predicted (see Fig. 4) while the 
AF2 phase corresponds to a spiral configuration in agree- 
ment with neutron scattering data [6]. According to our 
numerical calculation, the new intermediate phase AF3 
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TABLE I. Order parameters. 
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is associated with a collinear phase where only half of the 
moments order with S || b. As the field is applied, two 
spin-flop transitions (AFl-> HF1 and AF2 -> HF2) are 
found. The comparison of both phase diagrams indicates 
that the role of the anisotropic terms D z q and D xz q is to 
reduce the critical field of the AF1— >• HF1 transition, de- 
crease the stability range of the intermediate phase AF3, 
and lead to a new magnetic order HF3 in which half the 
moments align into the ac plane. These findings could ac- 
count for the fact that no spin-flop phase transition has 
been observed experimentally up to 16 T for the spiral 
phase AF2. 

Our principal conclusions are that a new collinear 
phase (AF3) has been detected by high resolution ultra- 
sonic velocity measurements which occurs between the 
paramagnetic and the previously identified spiral phase. 
The magnetic-field vs temperature phase diagram for B 
|| b has also been determined, revealing the existence of 
a new spin- flop phase (HF1). Complementary dielectric 
measurements also confirm that magnetoelectric effects 
only exist in the non-collinear phase. Verification that 
the new AF3 phase must exist is achieved by a Landau- 
type model based on rigorous symmetry arguments. Fur- 
thermore, the occurrence of such a collinear state, just 
above a non-collinear state, is confirmed in well studied 
frustrated RMn0 3 and RMn 2 5 systems [20-22], and the 
kagome compound M3V2O8 [23]. Finally, the proposed 
model accounts for the experimental phase diagram of 
CuO determined in this work and is potentially useful 
for the description of other monoclinic multiferroic sys- 
tems, in particular MnW0 4 [24] and AMSi 2 6 [25]. 



ACKNOWLEDGMENTS 

This work was supported by the Natural Science and 
Engineering Research Council of Canada (NSERC). 



FIG. 5. Magnetic field - temperature phase diagram of CuO 
for H || b derived from the Landau free energy. Dotted lines 
represent prediction with only one anisotropic term included, 
D y Q. The solid line is for the case where all anisotropic terms 
considered. 
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